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Abstract- Diesel engines are used in automotive and stationary 
applications. The main problem with diesel engines is emissions 
of nitrogen oxides (NOx) and particulates. In order to minimize 
the emissions, it is necessary to design the diesel engine with 
better in-cylinder flow (air-fuel mixing) and combustion 
process. Computational Fluid Dynamics (CFD) simulation 
helps to understand the Diesel engine temperature distribution 
and NOx species concentrations with respect to time. A small 
direct injection (DI) engine was chosen for the study. CFD 
simulation results were compared with that of engine emission 
tests. Results were found to be in agreement with NOx 
emissions. This paper also presents the simulation results of 
direct injection diesel engine in-cylinder flow (air-fuel mixing) 
and combustion. 
Index Terms - CFD, Diesel engine, combustion modeling 

I. INTRODUCTION 

Direct - injection diesel engines are used both in heavy- 
duty applications and light duty ones due to their high thermal 
efficiency and low CO, emissions [1]. Fuel economy and 
emission norms are compelling to develop more efficient 
and cleaner engines. It is necessary to improve the 
combustion process to reduce exhaust emissions. 
Computational Fluid Dynamics is popularly used in different 
stages of engine design and optimization. The combustion 
system performance can be better understood using CFD 
simulation. The challenge in using CFD is the complexity 
of interaction of flow, turbulence, spray and combustion 
inside the IC engine cylinder. High-pressure injection [2-5] 
and modification of combustion chamber geometry [6-9] 
were tried to reduce particulate emissions. Exhaust gas 
recirculation (EGR) [10] is another technique used to reduce 
NOx emissions. Supercharging coupled with better injection 
timing reduces PM emission and NOx reduction [11]. These 
methods however, have a trade-off relationship between NOx 
and particulate emissions. To optimize the combustion and 
emissions in diesel engine it is necessary to understand the 
flow inside the cylinder. This paper presents the CFD 
simulation of a diesel engine combustion process to predict 
the temperature distribution inside the cylinder and NOx 
emission with reference to crank angle variation in the 
expansion stroke. 

II. ENGINE PARAMETERS 

The specifications of the engine are given in Table I. The 
physical modeling of the injector was not used since preset 
models of fuel injectors were available in the software 
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TABLE I. 
ENGINE SPECIFICATIONS 



Model No 


4360 


Bore (mm) 


B2 


Stroke (mm) 


6B 


Displacement (cm J ) 


359 


Comprs sion ratiD 


1S:1 


EPM 


3600 


HP 


7.5 



used for analysis. The motion of valves also cannot be 
included in the simulation. 

III. COMPUTATIONAL DOMAIN 

The 3-Dimensional model shown in figure 1 is used for 
the analysis. The simulated engine has a direct injected diesel 
one with 4 - hole injector and eccentric bowl combustion 
chamber. For numerical integration the structured 
(hexahedral) mesh is used. During the solution, as the piston 
moves, the internal mesh structure deforms automatically to 
minimise the distortion of each individual cell. Figure 1 and 
3 shows the geometry and mesh distribution when piston at 
a bottom dead center (BDC) and the injection parameters 
are given in Table II. 






Piston ffi TDC 



-Computational Domain 



Figure 1 . Computational model of the Combustion chamber. 

The starting Crank angle when piston @ BDC is 180°. For 
each stroke corresponding crank angle is 1 80°, total for 4 
stroke 720° CA). Injection starts @ 343° CA means, when 
piston moves from BDC to TDC, the fuel injection starts 17° 
CA before TDC. Total injection duration 24° CA means fuel 
injection starts 17° CA before TDC and ends with 7° after 
TDC). Figure 2 shows the piston movement profiles. 

^AMAE 



AMAE Int. J. on Production and Industrial Engineering, Vol. 01, No. 01, Dec 2010 





Pistol promt 




= : 

E » 

E 

| « 

10 


i i i i 
^ Piston @ BDC 








: 






/ 








/ * 






/ 




\ 




/ 


^ 


\ 




/ 




\ 




~z 




\ 




/ 






r 


X 




\ 




/ 






V 


/_-™ ' 








\ y riaiun .U; iljl. 




1! 


B 20 £0 450 SO S30 7-50 S50 H 
Crank Ang b (D44O 



Figure 2. Piston movement profiles. 
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Section B 



Section A 
Figure 3. Computational mesh of the Combustion chamber. 

TABLE II. 
INJECTION PARAMETERS 



Typ- 


Direct inj ection with sinek 
central injector 


Numb =r f nozzles 


4 


TdIhI flow rate 
(e/cycle) 


0.05 


Start time (CA)0 


343 


Stop time (CA)0 


367 


Mean diameter (m) 


I e-06 



IV. METHOD USED FOR ANALYSIS OF 
SIMULATION DATA 

The numerical simulations in this investigation are 
performed based on the commercial computational fluid 
dynamics software. The various models used in the analysis 
to represent complex processes are explained below. 

A. Turbulent In-cy Under Flow Modeling 

In a DI Diesel engine, the fuel is sprayed directly into 
the cylinder through the fuel injection nozzle. The fuel is 
then broken into a number of droplets. These droplets 
undergo collision and coalescence processes, exchanging 
momentum and energy with the high temperature and 
pressure surrounding gases inside the cylinder. Finally the 
droplets evaporate into vapor and mix with air. The equations 
governing the in-cylinder flow are species transport 
equations, continuity equation, momentum equation, energy 
equation, gas state equation and turbulence equations. 
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Species Transport Equation: 
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source terms due to the chemical reaction and fuel spray, 
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d P rr ( 



- V 



c: 



\ 



it 



-.-) 



,■ 



Momentum Equation: 
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(3) 
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where ^ is the source term due to the fuel spray. 
Energy Equation: 

* ■■. ) • KpJ 

The last two terms on the right hand side of Eqn.(4) are 
the source terms due to the chemical heat release and fuel 
spray interactions, respectively. 
State Equation: 



:-0 



RT 



f— 1 



(5) 



The turbulence effect is accounted by using the standard 

k-^ □ model [1]. In this model, Reynolds stresses are related 

to the mean flow using the Bousinesq hypothesis. The 

turbulent kinetic energy k Hand turbulent dissipation rate <; 
I rj are solved through the rollowing* ~ 



lg transport equations. 
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where the values of the constants are given as follows: C 
=1.44, C , =1.92, C , = -1.0, Pr =1.0, Pr =1.3. W s is the 
source term due to interaction of air and fuel spray, which is 
the negative of the rate at which the turbulent eddies are 
doing work in dispersing the spray droplets. 
As the spray is the dominant feature of the flow with in a 
diesel engine combustion chamber, accurately resolving the 
spray atomization is a prerequisite for diesel engine 
simulation. Two mechanisms used for breakup of spray are 
Kelvin-Helmholtz (K-H) inertia instability and Rayleigh- 
Taylor (R-T) accelerative instability mechanisms. It is 
assumed that the mechanism responsible for the break up of 
the spray are classical fluid dynamic instabilities which act 
at the interface between two fluids of different densities. 
For both of the instabilities, the time that the spray droplet 
breaks up is determined from the growth rate of the fastest 
growing wavelength predicted by the classical instabilities. 
For the K-H mechanism, the wavelength of the fastest grow- 
ing wave is given by Reitz [12]. The R-T break up mecha- 
nism [14] is employed to account for the effect of rapid 
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deceleration of the drops on the atomization process. The 
K-H and R-T mechanisms compete to break up the spray. 
More details of break up mechanisms are described in 
Ref.[12] and [14]. 

B. NOx Model 

Both Nitric Oxide and Nitrogen dioxide are formed 
during the combustion process in a diesel engine. Generally 
speaking, NO, concentration is negligible compared with 
NO concentration. It is usually considered to be a transient 
intermediate species existing only at flame conditions, and 
it is subsequently converted back into NO in the post-flame 
region. Thermal NO formation is the most predominant in a 
diesel engine among the different NO formations. The most 
widely accepted simplified thermal NO formation 
mechanism is the extended Zeldowich-NO [19]. The 
elementary reactions for this mechanism are: 

:s;; 

(9) 
(10) 

Among the three reactions, the reaction (8) has the highest 
activation energy, due to the strong bond in the N^-molecule. 
It is the rate-limiting step of the thermal NO formation, and 
determines the overall process of the NO formation. The 
reaction rate coefficients for these reactions are given in 
Arrhenius form: 

{ID 



O-N; 


ZZNO + N 


N + Oj 


DDNO+O 


N + OH 


DDNO-H 



K , = AT 



sxp 



at ' /( mol — s ) 



The constants A, ? and E used for this investigation are from 
Hanson and Salimian, et al., [16]. They are given in Table 
III. 

TABLE III. 
CONSTANTS FOR THE KINETIC REACTION COEFFICIENTS. 



Equation 
Nd 

(B) 



Forward Reaction 



Backward Reaction 



A 
3x10' 



E 



A 



E 
:} 



V. RESULTS AND DISCUSSION 

The engine emission readings measured using gas 
analyzer and values are given in Table IV. 

TABLE IV 

GAS ANALYZER READING (AT FULL LOAD CONDITION AND AT 

3400 RPM) 



CO 

HC 

hexaie 

COj 



NO [ppm vjiU B43 

Temperature Distribution at various crank angles: 



[% keU 


3.71 


[ppmKfiU 


4539 


\% ssS 


7.(57 




C. Combustion Model 

The combustion inside the cylinder is modeled by four 
kinetic reactions and six equilibrium reactions. The four 
kinetic reactions include those three elementary reactions 
for thermal NO and a forward reaction representing the 
oxidation of the fuel, which is considered to be consisting 
of C p H, 6 only. The fuel oxidation reaction is given below. 
2 C 12 H^+ 370 2 ^24C0 2 + 26H 2 

D. Initial and Boundary conditions 

The computation started from the end of intake stroke at 
the bottom dead center corresponding to 1 80° Crank angle 
and terminated at the end of expansion stroke corresponding 
to 540° CA. The pressure and temperature of the air inside 
the chamber on the starting of the analysis were set as lOOkPa 
and 350K based on the values given by Bedford [17] and 
Strauss [18]. No slip conditions were used on the walls. 
Standard wall functions were used for turbulence models. 
The amount of fuel injected was 0.05 g/cycle. The normal 
operating parameters were : 
Injection duration: 24° CA; Injection start timing:343° CA 
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Fig 3. (e) CA - 375 deg. Fig 3. (f) CA- 385 deg. 

Figure 3(a) to 3(f) presents the in-cylinder temperature 
distributions from the CFD simulation of base case design. 
The temperature distributions plotted for different crank 
angles on a vertical plane. Figure 3(a) shows the maximum 
temperature of around 950K which is attained from 
compression stroke. From 343 - 720 crank angle the 
temperature distribution plotted in Figure 6. Figure 3 (c) 
shows the initial flame occurs just before compression stroke. 
The spray and flame reaches the edge of the piston bowl 
with in short period 7-10 CA deg. The concentration of 
oxygen in the high temperature zone is very high resulting 
in a high NOx rate of formation. Figure 4 and 5 shows the 
NO formation and Soot formation inside the cylinder. 



1 



' 



Figure- 4 . !MO formation. Fieurs 1«. Soot tormatiori. 

TABLE V 

[ppm KpVJ 
Experimait 843 

CFD Simulation 80 5 - B20 (max) 

The NO from emission test is 843 ppm at 3400 rpm. The 
detailed time history of spray, fuel mass fraction and 
temperature distributions provided by the CFD simulation 
are valuable towards gaining a better understanding of the 
features of combustion for given engine configurations. The 
temperature distribution inside the cylinder for different 
crank angles, and NO emissions are given above. NOx 
formation is highly sensitive to temperature and also effected 
by species concentration. In-fact the flame in the eccentric 
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bowl is not sufficient to burn a complete combustion because 
of the bowl shape. Table V presents the comparison of 
simulated and measured emission (NO) for the baseline case. 
The deviations from experiment and simulation results of 
NO emission are around 3-5 %. It was found that the general 
agreement between prediction and engine test was good. 
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CONCLUSION 

In-cylinder physical variations of a single cylinder diesel 
engine have been investigated using a CFD simulation code 
and validated with experimental measurements. Temperature 
profiles as well as NOx emissions contour from CFD has a 
good agreement with experimentally measured data. For 
reducing NOx emission it is no doubt that a precise estimation 
of flame temperature and a good description of reaction 
chemistry is essential. The detailed time history of spray, 
fuel mass fraction and temperature distributions provided 
by the CFD simulation is valuable towards gaining a better 
understanding of the feature of combustion. 

Since the adopted methodology for the base case has 
given a high confidence level for further optimization 
geometry (bowl shapes) and operating parameters (Injection 
duration, Start of Injection timings) can be carried out for 
reducing NOx emissions through CFD. 

NOMENCLATURE 

A = Constant 

E = Activation temperature 

e = Specific internal energy 

F = Rate of momentum gain per unit volume due 

to the fuel spray 

g = Gravity 

h = Specific enthalpy of species m 

q = Rate of heat release per unit volume 

q = Rate of heat gain per unit volume due to spray 

interactions 

R = Universal gas constant 

W = Molecular weight of species 

W = Source term due to the interaction of air and 

s 

fuel spray 
n = Mass density 
n = Mass density of species m 
n = Rate of mass gain per unit volume from 

chemical reactions 
n = Rate of mass gain per unit volume from fuel 

spray 
K = Reaction rate coefficient 



6 = Viscous stress tensor 

(j. = Viscosity 

i = Constant 

D = Diffusion CoefficientA. 
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